ON FOURIER TRANSFORMS OF RADIAL FUNCTIONS 
AND DISTRIBUTIONS 



LOUKAS GRAFAKOS AND GERALD TESCHL 

Abstract. We find a formula tliat relates the Fourier transform of a 
radial function on R" with the Fourier transform of the same function 
defined on R""'"^. This formula enables one to explicitly calculate the 
Fourier transform of any radial function f{r) in any dimension, provided 
one knows the Fourier transform of the one-dimensional function t H^> 
f{\t\) and the two-dimensional function {xi,X2) i— >■ f{\{xi, X2)\)- We 
prove analogous results for radial tempered distributions. 



1. Introduction 

The Fourier transform of an integrable function ^(x) on R" is defined by 
the convergent integral 

= I $(x)e-2'^"-«dx. 
Jr." 

If the function $ is radial, i.e., = for some function ip on the 

line, then its Fourier transform is also radial and we use the notation 

FnmO = M^)ir) , 

where r = |^|. In this article, we will show that there is a relationship be- 
tween J>^((/?)(r) and Tn+2i^)ii") as functions of the positive real variable r. 

We have the following result. 

Theorem 1.1. Let n > 1. Suppose that f is a function on the real line 
such that the functions f{\ ■ |) are in L^(Ti^^'^) and also in L^(R'^). Then 
we have 

(1) = 

Moreover, the following formula that is valid for all even Schwartz functions 
ip on the real line: 

(2) -Fn+2(^)(r) = ^^J'n{s'''^'^^{^{s)sn){r), r > 0. 
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2 L. GRAFAKOS AND G. TESCHL 

Using the fact that the Fourier transform is a unitary operator on L^(R") 
we may extend ([T]) to the case where the functions /(| -1) are in L^(R'^^^) and 
in L^(R"). Moreover, in Section[4]we extend ([T]) to tempered distributions. 
AppUcations are given in the last section. 

Corollary 1.2. Let f{r) he a function on [0, oo) and k some positive integer 
such the functions x — )■ are integrable over R" for all 1 <n < 2k + l. 

Then we have 

T (fV^ 1 ^ (-i)^(2fc-^-l)! 1 fdV 
^2k+i{f){p) - 1. 2^-^(/fc -£)!(£-!)! ^ [Tp J 

and 

T (fMn\- ^ V (-l)^(2fc -£-!)! 1 (dV 

The Corohary can be obtained using ([1]) by induction on k. The simple 
details are omitted. Again integrable can be replaced by square integrable. 

2. The proof 

The Fourier transform of an integrable radial function on R" is 

given by 

Mf)m=2^ f{s)(T-Y J^^l{2TTs\i\)sds 







oo 



=(27r)2 / f{s)Jn^{2TTs\i\)s^-'ds, 
Jo 

with u = ^ — 1 and Ju{x) = x~'^Jy{x). This formula can be found in many 
textbooks, and we refer to, e.g., O Sect. B.5] or [9J Sect. IV. 1] for a proof. 
Moreover, this formula makes sense for all integers n > 1, even n = 1, in 
which case 

, , [2 cost 

Let us set 

/•oo ^ 

n^-i{f){r) = (27r)t f{s)Jr,^^{2T,sr) s^-'ds . 

Then we make use of B.2.(l) in [3j, i.e. the identity 

(3) |:X(r) = -rX+i(r), 

which is also valid when v = —1/2, since 

' 2 sin t 
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In view of ([3]), it is straightforward to verify that 

provided / is such that interchanging differentiation with the integral defin- 
ing 'Hn_i is permissible. For this to happen, we need to have that 



\m\l{^-.{rs, 

and thus it will be sufficient to have 



^ds < oo 



(4) / \f{s)\rsV^{rs)\s^-'ds<c I 1/(^)1 " ' s^-'ds < oo 
Jo Jo {l + rs)~ 

since | J|(s)| < c(l + s)-"/2-i/2_ ^^^^^ . |) L^(W+^) we have 

fl/r POD 

n+1 ^ / . »/ M Il+i 
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(5) / / \f{s)\s—ds <oc 

Jo Jl/r 

and this certainly implies (j4j) for all r > 0. We conclude ([T]) whenever ([5]) 
holds. 

To prove this ([2j) we argue as follows. We have 

^n_i(r-+i-(^(rK))(r) = (2vr)t / Jn_,{2^sr)ds 

J 

and integrating by parts the preceding expression becomes 

/"OO ^ 

(27r)t+2 / (^(s)s"sr2 (27rsr)ds 
Jo 2 

which is equal to '2,'Kr^'H.r,A^_^{Lp){r^. This proves ([2]). 
Remark 2.1. A''oie that we have 

where 
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Hu{f){r) = fis)Mrs)sds, v> ^, 

is the Hankel transform. This of course ties in with the fact that the Hankel 
transform also arises naturally as the spectral transformation associated with 
the radial part of the Laplacian —A; we refer to [H Sect. 5] and the refer- 
ences therein for further information. Moreover, note that [6\ contains the 
associated recursion from Theorem for the Hankel transform, but only 
for even Schwartz functions. This recursion was rediscovered in connection 
with the radial Fourier transform in [8] for the case of Schwartz functions. 
See also [5\ for related results. 

A transference theorem for radial multipliers which exploits the connection 
between the Fourier transform of radial functions on R" and R"+2 ^Q^g 
obtained in [1]. This multiplier theorem is based on an identity dual to 



4 



L. GRAFAKOS AND G. TESCHL 



3. Radial distributions 

We denote by 5(R") the space of Schwartz functions on R" and by 5'(R") 
the space of tempered distributions on R". A Schwartz function is called 
radial if for all orthogonal transformations A G 0(n) (that is, for all rotations 
on R") we have 

(f = (f o A . 

We denote the set of all radial Schwartz functions by Srad{^")- For further 
background on radial distributions we refer to Treves [121 Lect. 5]. Observe 
that in the one-dimensional case the radial Schwartz functions are precisely 
the even Schwartz functions, that is: 

5rad(R) = SevenCR) = W ^ ^(R) : ip{x) = (p{-x)}. 

Similarly, a distribution u G 5'(R") is called radial if for all orthogonal 
transformations A G 0{n) we have 

u = u o A . 

This means that 

{u,ip) = {u,ip o A) 

for all Schwartz functions (/? on R". We denote by 5',^^(R") the space of all 
radial tempered distributions on R". We also denote by S"""^ the (n — 1)- 
dimensional unit sphere on R" and by oJn-i its surface area. 

Given a general, non necessarily radial, Schwartz function there is a nat- 
ural homomorphism 

5(R") ^ SradCR), ^(x) ^ (/p°(r) = [ ^{rO) dO 

with the understanding that when n = 1, then ip"{x) = ^((/^(a::) -|- (p{—x)). 
Conversely, given an even Schwartz function on R we can define a corre- 
sponding radial Schwartz function via 

5,arf(R) ^ 5,arf(R"), V9(r) ^ ^^{x) = V9(|x|). 

The map i— )• ip^ is a homomorphism; the proof of this fact is omitted since 
a stronger statement is proved at the end of this section. Both facts require 
the following lemma: 

Lemma 3.1. Suppose that f is a smooth even function on R. Then there 
is a smooth function g on the real line such that 

fix) = g{x^) 

for all X G R. Moreover, one has for t > 

(6) \g^'Ht)\<C{k) sup \f^^'\s)\. 

0<s<v^ 
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Proof. By Whitney's theorem \13\ . there is a smooth function g on the real 
hne such that 

fit) = 9{t') 

for all real t. 

To see the last assertion we use the following representation of the re- 
mainder in Taylor's theorem: 




which yields in particular: 

k\ {2k)\ 

Using these identities, one easily derives ([6]). □ 

The composition ip i— t- {ip°)'^ = (p'^"''^ gives rise to a homomorphism from 
5(R") — 5rad(R-") which reduces to the identity on radial Schwarz func- 
tions. In particular, the map cp ip° defines a one-to-one correspondence 
between radial Schwartz functions on and even Schwartz functions on 
the real line. Moreover, ip is radial if and only if 99 = tp^'"-'^. 

Proposition 3.2. For u G 5^^^(R") and p G 5(R") we have 

(U,(^) = (^Z,(^™'^). 

Proof. By a simple change of variables the formula holds for any u which is 
a polynomially bounded locally integrable function. Next we fix a tempered 
distribution u on R" and we consider a radial Schwartz function ip with 
integral 1 and we note that the convolution ol ips * u converges to u in 
5'(R"') as e — )• 0. Hence, since the claim holds if u is replaced by * u by 
the first observation, it remains true in the limit e — )• 0. □ 

In particular, note that a radial distribution is uniquely determined by 
its action on radial Schwartz functions. Furthermore, given a distribution 
u G 5'(R") we can define a radial distribution u^"''^ G S'^^^CR"') via 

{u'^'^'^^ip) := (u,v9™'^). 

Moreover, u is radial if and only if n = u^"''^. 

For n G Z+ we denote by TZn = r"~^5ei,en(R) the space of functions of 
the form ^(r)r"~^, where ip is an even Schwartz function on the line. This 
space inherits the topology of S'(R) and its dual space is denoted by TZ'^. 
Two distributions tfi , ^2 G S' (R) are equal in the space 7^'„ if for all even 
Schwartz functions ^ on the line we have: 

(wi,r"'"V('')) = (w2,r'"'"V('")) • 

Note that in dimension n > 2 we have that all distributions of order n — 2 
supported at the origin equal the zero distribution in the space TZ'^^. Thus 
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two radial distributions wi and 102 are equal in TZ'^ whenever wi — 102 is a 
sum of derivatives of the Dirac mass at the origin of order at most n — 2. 

One may build radial distributions on R*^ starting from distributions in 
TZ'^. Indeed, given in 7^'„ and ip in 5(R'^) we define a radial distribution 
u by setting 

{u,(p) := (uo,(/? (r)r ) 

The converse is the content of the following proposition. 

Proposition 3.3. The map IZn — >• Sradi^"'), ^{r)r^~^ 1— )■ 'il>'~'{x) is a home- 
omorphism and hence for every radial distribution u we can define in TZ'^ 
via 

(no,^(r)r"-i) := -^{u,i;''). 

Proof. It suffices to show the first claim. To this end we will show that for 
all multiindices a and (3 we have 

sup < E sup|r-(-^)'(r«-V(r-))|. 

tpR" — r>n \ar / 

0<i,m<i{\P\+\cx\+n)^^ 

First we consider the case |x| < 1. Setting r = |x| < 1 we have 

m m 

k=0 k=0 

m m 

k=0 k=0 °<"<^ 

using Lemma |3. 11 

We will make use of the inequality 

.It) 

which follows by applying the fundamental theorem of calculus M times and 
of the identity: 

(«) >"S^(«)^i:(-»'«(7)G0(|)"-(^"-'«»)) 

which is valid for M > m and is easily proved by induction. 
Applying ([7]) to ■ip^'^''\s) we obtain 

(9) |^(^^)(.)| < sup 

0<t<s ^atJ 

and using ^ for s^^V'(2^)(s) with M = 2|/3| + n - 1 and m = 2k we de- 
duce that \^^'^^\s)\ is pointwise bounded by a sum of derivatives of terms 
ip{^) multiplied by powers of s. It follows that sup^^^Q IS con- 

trolled by a finite sum of Schwartz seminorms of the function s^~^ip{s). 



(7) |^(.)| < sup (-) {t^'m){s 

0<t<s ^ot/ 
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The case \x\ > 1 is easier since when / 

m ^ 

and taking M = max(|a|, \ j3\ + n — 1) we have 

(10) sup <C^^sup{s^V^^'^(^)l}, 

I^I>1 ~t«>l 

which is certainly controhed by a finite sum of Schwartz seminorms of 
in view of □ 



Note that if u is given by a function /(x), then is given by the function 

l_d_ 

r dr 



f°{x). We also remark that the map -4- is a homomorphism from IZ'^ to 



Tl'^j^i defined as the dual map of — 

A related approach defining for a given distribution u supported in 
R" \ {0} can be found in [10]. Our approach does not impose restrictions 
on the support of the distribution. 

4. The extension to tempered distributions 

Let ti be a radial distribution on and let -Fjt(n) be the /c-dimensional 
Fourier transform of u. 

Theorem 4.1. Given an even tempered distribution vq on the real line, 
define radial distributions Vn on R" and Vn+2 on R""*"'^ via the identities 

(11) {vn, (f) = {vo, ia;„_ir"~ V) 

for all radial Schwartz functions f{x) = ip"{\x\) on R" and 

for all radial Schwartz functions <f{x) = ip°{\x\) on R""*"'^. 

Now consider n" = Fn{vn) o-nd = Fn+2{vn+2) ■ Then the identity 

(12) -:t^ 



27rr dr 



holds on ■ 



Proof. We denote by (•, •)n the action of the distribution on a function in 
dimension n. Let ip{r) be an even Schwartz function on the real line. Then 
we need to show that 

(13) ( - ^|:<,u;„+ir"+V(r))^ = u;„+ir"+V(r)>i . 
This is equivalent to showing that 

(14) i^«,..„+i(rXr))')^ = «+2^^„+ir"+V(r)>i . 
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We introduce the even Schwartz function r]{r) = r~'^^^{r^Tp{r)y = ntp^r) + 
rip'{r) on the real hne and functions t]'~^ on R" and ip'^ on R""'"^ by setting 

V'^(x) = v(|x|) v''{y) = vi\y\) 

for y e R" and x G R"+^. Then (fill) is equivalent to 

(15) ^^«''^"-i"""'^(^)>i = «+',^n+ir"+VW>, 
which is in turn equivalent to 

(16) ^Z^(^-^--)'^'')n = {FnMv^+2UO^^^, 

and also to 

(17) ^^^-^"(^'')>. = (-"+2'^"+2(^^)>„+2- 

We now switch to dimension one by writing ([T7]) equivalently as 

(18) -l^(^;o,a;„_ir"-ij-„(77)(r)>^ = (^;o, ^n+ir"+i J-„+2(V')(r)>, . 
But this identity holds if 

l-Tn{r,){r)=r^Tn+2mr). 



which is valid as a restatement of ([2]); recall that T/(r) = r "'*'"^^(r"^/;(r)). 
This proves ([13]). □ 

It is straightforward to check that for polynomially bounded smooth func- 
tions all operations coincide with the usual ones. We end this section with 
a few more illustrative examples. Let (5„ be the Dirac mass on R". 
Examples: 

a) Let Vn = Sn- One can see that 



^° ujn-i{n - 1)1 ydr J 



satisfies (fTTI) . Acting vq on r"+^(^°(r) yields that Vn+2 = and thus uj"'"^ 
0. Also u" = 1; so both sides of (fT2l) are equal to zero, 
b) Let Vn+2 = Sn+2- Then 



ujn+i{n + 1)! \dr/ 



Let A. = df + ■ ■ ■ + he the Laplacian. We claim that the distribution 

(19) Vn = ^-A{5n) 

UJn+1 n 

satisfies ([H]). Then = 1 and also n" = -r^(27r)^a;„_i/(2ntj„+i). Thus 
(fT2]) is valid since 2iTUJn-i = nuin+i- 
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It remains to prove that the distribution Vn in ()19p satisfies (jlip . For 
(p{x) = ip°{\x\) in 5(R") we have 

(20) {vn,v) = (^o,a;„-ir"-V°(r)> = "^"'^ ^ , „ (^i, V°(0)^""^^) 

w„+i(n + l)! 

and one notices that the (n — l)st derivative of r"'~^ip°{r) evaluated at zero 
is equal to ^{n + l)l{ip°)"{0). To compute the value of this derivative we use 
Lemma O to write ip{x) = ip°i\x\) = g{\x\'^) where g'{0) = i(y?°)"(0). It 
follows that ^'(0) = 2^A((/3)(0). Combining these observations yields that 
the expression in ([20]) is equal to 

^-aM(o) = (^^1a(o,9.), 

LOn+in \ujn+in I 

which proves the claim. 

Remark 4.2. As 'pointed out in Remark \2.1{ the action of the Fourier trans- 
form on the associated function on the reals 99° is given by the Hankel trans- 
form. In particular, the results in this section also give a natural extension 
of the Hankel transform (for half-integer order) to distributions. Of course 
this coincides with the usual approach, see [H [HI [T5] and the references 
therein. To this end observe that the space F used in [6] is precisely the set 
of functions on [0,oo) which extend to an even Schwartz function on R. 

5. Applications 

We begin with a simple example. In dimension one we have that the 
Fourier transform of sech(7r|x|) is sech(7r|^|). It follows from ([T|) that in 
dimension three we have 

F3(sech(^|x|))(C) = ^sech(7r|e|)tanh(7r|e|). 

since 

d 2 e'^^ — e~'^^ 1 

— = —2tt- -TT = — 27r— sechfvrr) tanhfvrr) 

dr e^' + e-'^^' (e'^'' + e-'^'')^ 2 ^ ' ^ ' 

Continuing this process, one can explicitly calculate the Fourier transform 
of sech(7r|x|) in all odd dimensions. 

More sophisticated applications of our formulas appear in computations 
of functions of the Laplacian —A, which arise in numerous applications. 
For example, in quantum mechanics the Laplacian —A arises as the free 
Schrodinger operator (cf., e.g., [7], pT]) and functions /(—A) are defined 
via the spectral theorem by 

/(-A)^ = K*^, ¥'G5(R"), 

where K is the tempered distribution given by the inverse Fourier transform 
of the radial function /(47r^|^p), which is assumed polynomially bounded. 
Knowledge of the inverse Fourier transform of /(47r^|^p), for ^ G R and 
^ G R^, yields explicit formulas for the kernel K of /(—A) in all dimensions. 
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An important application is the explicit calculation of the n-dimensional 
kernel Gn{x) for the resolvent associated with the function /(r) = {r — z)~^, 
z € C\[0,oo). In the one-dimensional case, an easy computation shows that 

= ^^e-^l^l. 



Hence, by the I? version of Theorem ll.il (cf. the discussion right after The- 
orem [TTT]) the three-dimensional kernel is given by 



-z\x\ 



--\x\ 47r|a;| 

The computation of ^5(2;), G7(x), . . . requires Theorem 14.11 since the as- 
sumptions of Theorem 1 1.1 1 are no longer satisfied. For instance. Theorem 14. II 
gives 

Another interesting situation where our theorem is useful are the spectral 
projections associated with the function /(r) = X[o,s]('^)) -^^ > 0. Again in 
the one-dimensional case the kernel for the resolvent can be easily computed 
and found to be 

TTX 

Thus by Theorem II. II the three-dimensional kernel is given by 

Id,, sm(\x\^/E) — |x|\/^cos(|x|\/^) 

-^3(2;) = -- — , , = o 21 13 • 

Zirr ar r=\x\ 2tt'^\x\-^ 

Finally, the Fourier transform is a crucial tool in solving constant coef- 
ficient linear partial differential equations (cf., e.g, [2]). Using the above 
trick one can of course derive the fundamental solution for the heat (or 
Schrodinger) equation in three dimensions from the one-dimensional one. 
However, since the three-dimensional case is no more difficult than the one- 
dimensional case we rather turn to the Cauchy problem for the wave equation 

utt - Au = 0, n(0, x) = ipix), ut{0, x) = v?(x), 

in R", whose solution is given by 



u{t,x) = cos{t^/^)^/J{x) + ip{x). 

V-A 

Since the first term can be obtained by differentiating the second (with 
respect to t) it suffices to look only at the second and assume -0 = 0. 
Moreover, since the Fourier transform of /(x) = is Fi{f){^) = 

|ar^X[-i/2,i/2](C/«), we obtain 

u{t,x)= I ^X[-t,t]{x -y)v{y)dy, 
Jn ^ 
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which is of course just d'Alembert's formula. In order to apply Theorem 14. II 
we use vo{r) = such that = F^^{vi) as well as ul are associated 

with the function ^X[-t,t](.x)- Hence by Theorem 14.11 

and we obtain Kirchhoff's formula 

u{t,x) = ^ (p{x- te)de. 
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